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Interest in non-standard recombination scenarios has been spurred by recent cosmic microwave 
background (CMB) results from BOOMERANG and MAXIMA, which show an unexpectedly low 
second acoustic peak, resulting in a best-fit baryon density that is 50% larger than the prediction 
of big-bang nucleosynthesis (BBN). This apparent discrepancy can be avoided if the universe has 
a non-standard ionization history in which the recombination of hydrogen is significantly delayed 
relative to the standard model. While future CMB observations may eliminate this discrepancy, 
it is useful to develop a general framework for analyzing non-standard ionization histories. We 
develop such a framework, examining non-standard models in which the hydrogen binding energy 
Eb and the overall expression for the time rate of change of the ionized fraction of electrons are 
multiplied by arbitrary factors. This set of models includes a number of previously-proposed models 
as special cases. We find a wide range of models with delayed recombination that are able to fit 
the CMB data with a baryon density in accordance with BBN, but there are even allowed models 
with earlier recombination than in the standard model. A generic prediction of these models is 
that the third acoustic CMB peak should be very low relative to what is found in the standard 
model. This is the case even for the models with earlier recombination than in the standard model, 
because here the third peak is lowered by an increased diffusion damping at recombination relative 
to the standard model. Interestingly, the specific height of the third peak depends sensitively on the 
model parameters, so that future CMB measurements will be able to distinguish between different 
non-standard recombination scenarios. 



I. INTRODUCTION 

In the past year, observations of the cosmic microwave background (CMB) fluctuations by the BOOMERANG § 
and MAXIMA |^] experiments have produced data of unprecedented precision on CMB fluctuations at smaU angular 
scales. While generally confirming the adiabatic, flat {fl « 1) model predicted by inflation, these observations have 
several puzzling features. In particular, the position of the first acoustic peak is at a slightly larger angular scale than 
is predicted in the fiat model, and the amplitude of the second peak is unexpectedly low (see, e.g., Ref. |^). If these 
results are fit using the standard set of cosmological parameters, the result is a CMB prediction for the baryon density 
of ilft/i-^ ^ 0.03 Q - 1^. In contrast, the prediction for ilf,/i^ from Big-Bang nucleosynthesis is 51f,/i^ ^ 0.02 ||^,||. 

This apparent discrepancy could easily vanish in the light of future CMB measurements. For the time being, 
however, it has led to a great deal of interest in models with non-standard ionization histories, since one way to 
explain the CMB observations and preserve agreement with the BBN baryon density is to postulate that the epoch 
of recombination was delayed to a lower redshift than in the standard model. Peebles, Seager, and Hu suggested that 
this could occur if there were sources of Ly a photons present at the epoch of recombination p\. A more speculative 
mechanism is a time- variation in the fine-structure constant, a pO|- p^ . The authors of Ref. 1^ used such a time- 
variation (along with changes in the cosmological parameters) as an example of how one might model non-standard 
recombination in general. 

Because there are a variety of models with non-standard ionization histories, we feel that is it worthwhile to try to 
develop a general framework for analyzing such models. Ideally, one would like to investigate the consequences for the 
CMB of an arbitary Xe{z), the ionization fraction as a function of redshift. It is obviously impractical to investigate 
arbitrary functional forms for Xe- Instead, we have attempted to parametrize deviations of Xe{z) from the standard 
ionization history in terms of a small number of physically-motivated parameters. In particular, we multiply the overall 
ionization/recombination rates by a free parameter a, and the binding energies by a parameter b. We then calculate 
the predictions for the CMB fiuctuation spectrum as a function of these parameters and compare with observations. 
Our parametrization is discussed in the next section, along with the results of modifying the ionization history. Our 
conclusions are given in Sec. 3. Although we have compared our results with the recent BOOMERANG and MAXIMA 
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experiments, our work provides a general framework for discussing non-standard recombination scenarios, and it can 
be applied to future observations as well. 



II. NON-STANDARD RECOMBINATION 



Given that one cannot examine all possible recombination histories, what is the best subset to investigate? We 
have attempted to model a subset of such recombination histories which is physically motivated and, at least to some 
extent, reduces to the Peebles et al. model Q and the time-varying a model as special cases. There are 

two possible approaches to modifying the ionization history Xe{z). We can directly modify this ionization history, 
or we can alter the evolution equation for dxe/dt. We use the latter approach, but we also examine the relationship 
between variations in dxe/dt and the resulting form for Xe{z). 

Consider first the equation for dxe/dt: ||l7|] : 



dxe 
~dt 



TZupxl — (3{1 — Xe) exp 



Bi — B2 
kT 



(1) 



where TZ is the recombination coefficient, (3 is the ionization coefficient, Bn is the binding energy of the n'^ H-atom 
level and Up is the sum of free protons and H-atoms. The Peebles correction factor (C) accounts for the effect of the 
presence of non-thermal Lyman-a resonance photons; it is defined as 



C 



1 + KA{1 



1+K{A + I3){l-Xe)' 



(2) 



In the above, K = H^^ripC^ / 8^1/12 (where 1^12 is the Lyman-a transition frequency), and A is the rate of decay of the 
2s excited state to the ground state via 2 photons The ratio (3/TZ is fixed by detailed balance. 

We modify the evolution history for x^ as follows. We introduce two new parameters, a, and b, into equation (Q). 
The parameter a multiplies the overall rate for dx^/dt, while b multiplies all of the binding energies i3„. Then equation 
(lil) becomes: 



dXe 

~dt 



aC 



TZnpx\ — (3{1 — Xe) exp —6 



B i — B2 
kT 



(3) 



We take a and b to be constants independent of z (a model with redshift-dependent a and b would simply take us back 
to arbitrary functional behavior for Xe{z)). The case a = b = 1 corresponds to the standard model. In order to be 
completely self-consistent the recombination rate for helium should also be changed with varying a and b. However, 
this is a very small effect and in all our calculations the recombination history of helium is assumed to follow the 
standard model. This assumption has no bearing on any of our conclusions. 

Our a and b parameters have a simple physical interpretation. A change in a alone represents a change in the 
ionization and recombination rates which preserves detailed balance at a fixed binding energy, since the ratio of the 
ionization to recombination rate is unchanged. A change in b alone simply shifts the epoch of recombination up or 
down in redshift by a fixed amount. 

This particular parametrization has several advantages. The time- varying a model is basically a special case of this 
model 11|. In the Peebles et al. model |^ delayed recombination arises from additional Ly-a resonance photons 
produced by some unknown source. We have not directly tested this specific model in (a,5)-space, but instead we have 
examined a related model in which the rates for the Ly-a and 2s ^ Is transitions to the ground state are reduced by 
a factor e which is assumed to be constant in time. In this model the Peebles correction factor is changed relative to 
Eq. (0) 



C 



l + KK{l 



l + K{K + p/e){l-Xe)' 



(4) 



Although this model is very similar in behaviour to that of the Peebles et al. model Q, they are not completely 
equivalent. However, the exact form of the Peebles et al. model is quite speculative, and our purpose here has just 
been to test some of the possible ways to alter recombination. The Peebles et al. model could be mimicked almost 
exactly by a time-dependent e, but this is an unnecessary complication. 

In order to examine the viable region in (a, b) parameter space we have performed a analysis on the data from 
the BOOMERANG 0] and MAXIMA [| experiments. Our procedure is to maximize the likelihood for each point in 



2 



(a, 6) space for the following free parameters: the total matter density, fim, the Hubble parameter, Hq, the spectral 
index of the primordial power spectrum, rig, and the overall normalization of the spectrum, Q. We have assumed a 
flat geometry so that fl\ = 1 — ilm, and also that reionization is unimportant (assuming that the optical depth to 
reionization is small, t ~ 0). The vector of free parameters is then 

9^{a^,Ho,n„Q}, (5) 

and the likelihood function to be optimized is 



\2 



£cxAexp( ). (6) 





FIG. 1. The allowed region in a,b parameter space from the combined BOOMERANG and MAXIMA data, assuming a 
baryon density of Q,t,h^ = 0.019. The dark shaded (green) region is the la allowed region, and the light shaded (yellow) 
is the 2a region. Also shown: la) The curve in a, b space which best fits the varying fine-structure constant (a) model for 
non-standard recombination; Aa is defined as Aa = (a — ao)/ao, lb) The curve which best fits the modified Peebles et al. 
model for delayed recombination. The quantity e is defined in Eq. (^. 

The agreement between our models and the observations, as a function of a and 6, is shown in Figs, la and lb, 
for ilfc/i^ = 0.019, and in Fig. 2 for flth^ — 0.03. In both figures, the area outside the dark (green) shaded region is 



3 



excluded at the 1 — a level, and the area outside the light (yellow) shaded region is excluded at the 2 — a level. The 
allowed region is a broad band, in which larger values of a are compensated by smaller values of b. (Note that the 
apparent cutoffs at small a and small b in Fig. 1 are the real boundaries of the confidence region) . This result hides 
a great deal of information: from these graphs alone it is impossible to tell whether the allowed region corresponds 
to identical forms for Xe{z) produced by different values of a, b, or whether Xe{z) varies greatly within the allowed 
region. The latter is, in fact, the case. This is shown in Fig. 3, in which we graph the ionization history for three 
pairs of a, b within the allowed region. Increasing b and decreasing a results in a surface of last scattering at higher 
redshift which is much broader. 

We also display, for comparison, the values of a and b which best fit the models with a change in a (Fig. la) and 
our version of the Peebles et al. model, given by equation (^ (Fig- lb). The behavior of the time- varying a model is 
easy to understand; an increase in a results in an increase in all of the binding energies (and thus, an increase in b) 
and it also increases the ionization and recombination rates (an increase in a) |10 1^. Thus, the curve corresponding 
to time-varying a runs almost perpendicular to our best-fit contour, and we find good agreement for negative values 
of Aa, as in Refs. - In Fig. lb, we see that the correspondence between the modified Peebles, et al., model 
and our a, b formalism is more complicated. The part of the curve corresponding to e < 1 corresponds to faster than 
normal recombination and is not physically related to the Peebles et al. model. For e > 1 both a and b are changing 
with changing e. From Fig. 1 in Ref. it is possible to understand the path taken by the curve in a, 6-space. As 
e is increased recombination is pushed to smaller z, but at the same time the width of the recombination surface 
decreases. Therefore, for increasing e, the best-fit curve should move down and to the right in a, b space, exactly 
what is seen from Fig. lb. This is different from the varying a curve, where the width of the recombination surface 
increases as a decreases. This means that the best-fit values in the time-varying a model and the best-fit values in 
the modified Peebles model lie in slightly different parts of the a, b plane. 
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FIG. 2. The light (yellow) and dark (green) shaded regions are equivalent to those in Fig. 1, except that they are calculated 
assuming a high baryon density, Q,bh? = 0.030. 
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FIG. 3. Three different ionization histories, aU chosen to lie within the aUowed region. The dotted curve is for a — 0.1, b — 1.0, 
the sohd curve for a — 1.0, b — 0.75 and the dashed for a = 10.0, b = 0.65. 

On first sight the results in Fig. 3 would mean that the large allowed region in (a, 6)-space is entirely due to 
degeneracy between a, b and the other cosmological parameters. In order to investigate this possibility we have 
performed a simple Fisher matrix analysis of how degenerate a and b are with the other cosmological parameters. 
The Fisher information matrix is given by 



where i and j denote elements in the vector of cosmological parameters to be determined and cr{Ci) is the uncertainty 
in the measurement of Ci. The standard deviation in a measurement of parameter Oi is then given by a{0i)'^ = {F~^)ii. 
For simplicity we assume that the measurement error in the C; 's is purely due to cosmic variance so that 



The outcome of this analysis is shown in Figs. 4a and 4b. The precision with which either a or & can be determined 
drops considerably if a and b must be determined simultaneously, meaning that a and b are to some extent degenerate. 
However, the biggest loss of precision still comes when the other cosmological parameters have to be determined as 
well. The partial degeneracy between a and b has a simple physical explanation. If the redshift of recombination is 
lowered the radiation content at recombination is lower, leading to a smaller early ISW effect. On the other hand, if 
the width of the recombination surface is narrowed, the diffusion damping of fluctuations is smaller. These two effects 
can to some extent compensate each other, and lead to a partial degeneracy between a and b. However, the degeneracy 
is not exact, because the early ISW effect dominates near the first peak in the spectrum, while the diffusion damping 
dominates at larger I. 
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FIG. 4. The expected accuracy with which a can be measured, as a function of the maximum I that the power spectrum can 
be measured to. We have assumed that the power spectrum can be measured at all /, up to imax, to within cosmic variance 
a{Ci)/Ci = [2/{2l + 1)]^/^ Fig. 4a) The sohd hue is the expected precision if a alone needs to be determined from the 
CMB, the dotted is the case where a and h must be simultaneously determined, and the dashed the case where all parameters 
{Qm,^b, HojTisjQ, a,b) must be determined. Fig. 4b) The same as for Fig. 4a, except that a and 6 are interchanged. The 
fiducial model is a = 6 = 1. 



A second possible approach to the problem of the CMB as a function of a generic recombination history would be 
to modify Xe{z) directly by hand. This is a somewhat unphysical approach, so we have instead examined the general 
behavior of Xe{z) as a function of a and b. The two key properties of the ionization history which affect the CMB 
are the redshift of last-scattering, and the width of the last-scattering surface. We parametrize the former in terms 
of ^1/2, the redshift at which Xe = 0.5. To estimate the width of the last scattering surface, we define the parameter 
Az to be z{xe = 0.9) — z{xe = 0.1). While both of these definitions are somewhat arbitrary, they serve the desired 
purpose of indicating the redshift and width of the last scattering surface. The behavior of these quantities as a 
function of a and b is illustrated in Figs. 5 and 6, and in Fig. 7 we show Az/zi/2 as a function of a and b. 
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FIG. 5. The redshift of last scattering (21/2) as a function of a and h. 
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FIG. 6. The width of the last scattering surface (Az) as a function of a and h. 
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FIG. 7. The width of the last scattering surface relative to the redshift of last scattering, Az/ zi/2, as a function of a and 6. 

Several features are obvious from these plots. In Fig. 5, we see that Zi/2 is a function only of h (and is independent 
of a) in the limit of large a. We see the same effect in Fig. 6, in which Az becomes a function only of 6, but only for 
the case of small b. These results follow from the fact that for large a, the ionization fraction tracks its equilibrium 
value nearly exactly, and the equilibrium value of Xe at fixed z depends only on the binding energies (and hence, on 
h). This argument breaks down for small a, since in this case the ionization fraction no longer tracks the equilibrium 
abundance. The argument breaks down for large h for a more complicated reason. Roughly speaking, Xe tracks its 
equilibrium abundance as long as (dxe/dt)/ H > 1, where H is the expansion rate. At large b, recombination occurs 
earlier, when H is larger (smaller cosmic time), requiring a larger value of a to maintain equilibrium. This effect is 
more apparent in Fig. 6 than in Fig. 5 because of the way that we have defined Zi/2 and Az. The value of Zi/2 is 
the redshift at which Xe = 0.5, while Az depends on the much later redshift at which Xe — 0.1. Hence, zi/2 will be 
independent of a as long as equilibrium is maintained down to the redshift at which x^ — 0.5, while in order for Az 
to be independent of a, equilibrium must be maintained down to the much lower redshift at which Xe = 0.1. 

In the limit of extremely small a, recombination is continuing at the present, and the lower limit x^ = 0.1 we use 
to calculate Az is never reached. Hence, there is a limiting value for a below which Az is undefined; this is reflected 
in the degeneracy of our contours for small a in Fig. 6. 

Finally, if we define the width of our last scattering surface relative to the redshift of last scattering (as in Fig. 
7), then we see that Az/zi/2 is essentially a function of a, and is nearly independent of b. At large a, the ionization 
fraction tracks its equilibrium value nearly exactly (and so is independent of a). Furthermore, when Xe{z) is given by 
the Saha equation, Az/zi/2 is independent of the binding energy. Thus Az/zx/2 becomes independent of both a and 
b for large values of a, as seen in Fig. 7. 

Figs. 5 and 6 indicate that scanning over all possible values of a and b is nearly equivalent to scanning over all 
possible values of the redshift and width of the last scattering surface. Hence, our method of varying a and b provides 
a quite general study of arbitrary ionization histories. The one exception is the case of large zi/2 and small Az, which 
does not correspond to any values of a and b. However, this case corresponds to a decrease in Xe which occurs faster 
than for the equilibrium case. It is difficult to imagine a mechanism for achieving this. Furthermore, we expect the 
Ci spectrum to become independent of Az for sufficiently small Az, since the thickness of the last-scattering surface 
will become irrelevant when it becomes so narrow that diffusion damping can be neglected. For completeness we have 
also calculated for the case of Az = 0, for which we take Xg to be a step function. These are shown in Fig. 8. We 
see that no good fit is obtained for any value of Zi/2 (Xmin/'^-O-f — 3.8). 

If we consider the allowed region in Fig. 1 as a function of Zi/2 and Az, rather than a, and 6, we find that both Az 
and zi/2 are restricted by the current CMB observations. At the 95% confidence level, we have 300 < Az < 900 and 
500 < zi/2 < 2000. From Fig. 7, we see that our limits translate into an upper bound on Az/zi/2 of Az/zi/2 < 0.5. 
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FIG. 8. The x for the case of Az = as a function of zi/2- 




FIG. 9. The height of the first peak relative to the height of the third peak, X = Ci( 1st peak)/C; (3rd peak), as a function 
of a and b. The dashed curve gives the standard model {a = b = 1) value oi X = 2.65. 



III. CONCLUSIONS 



Our investigation of the effect of a non-standard ionization history on the CMB using a change in the overall 
ionization/recombination rates (through the parameter a) and in the binding energies (through the parameter h) 
appears to provide a very general framework for studying such variations in the standard model. In particular, it 
seems possible to model arbitrary changes in both the epoch of last scattering and in the width of the last-scattering 
surface through such variations, with the exception of ionization histories having a narrow, high-redshift surface of 
last scattering, for which the ionization fraction decreases faster than in the equilibrium case. 



9 



Our results indicate that the BOOMERANG and MAXIMA results with n^h'^ = 0.019 (from BBN) can be well-fit 
over a broad range of choices for a and b. However, the lower bounds on a and b are quite robust (i.e., nearly 
independent of variation in the other parameter). We find that b > 0.5 and log(a) > —1.6 at the 95% confidence level. 
If we restrict a to its standard model value, changing only the binding energies, the allowed (95%) region is given 
by 0.5 < b < 1.1. Similarly, if we restrict the binding energies to be unchanged {b = 1) and change only the overall 
rates, then a is constrained to lie in the range — 1.6 < log(a) < 0.4. Our best- fit model for ^l^h^ — 0.019 has a per 
degree of freedom of 0.50. In comparison, the best fit for standard reionization {a — b — 1) with flbh^ — 0.030 gives 
per degree of freedom of 0.82, so both models are good fits. (These numbers include a 10% calibration error for 
BOOMERANG and 4% for MAXIMA). 

One robust prediction of all of these models is a decrease in the amplitude of the third peak, relative to its height 
in the standard model. In Fig. 9, we show the ratio X of the height of the first peak to the height of the third peak 
as a function of a and b. The standard model {a = b = 1, fl^h^ — 0.019) value, X — 2.65, is shown as a dashed 
curve. The entire la region gives values for X larger than in the standard model. In comparison, the best fit model 
with high baryon density {Vli,h? = 0.030) and a standard recombination history produces a much higher third peak 
{X = 1.82) than in the standard model. 

Our results are more general than previous comments that models with late recombination can be distinguished 
from high Vl}jh^ models by the amplitude of the third peak. Our allowed region includes models with recombination 
at slightly higher redshift than in the standard model (these models lie on the extreme left-hand side of Fig. 9). For 
these models, the third peak is still reduced in amplitude, but this reduction is due to diffusion damping from an 
increase in the width of the surface of last scattering, rather than from a decrease in the redshift of last scattering 



as in Refs. - 15 1. Our results seem to imply that any modification in the recombination history which fits the 
BOOMERANG and MAXIMA observations will result generically in a decrease in the height of the third peak. 

However, we also find that the actual height of the third peak is a function of a and b. From Fig. 9 it can be seen 
that this ratio increases as we move through the allowed region from small a, large b to large a, small b. So if the 
universe did have a non-standard ionization history, the amplitude of the third peak as shown in Fig. 9 should allow 
us to determine precisely the nature of the deviation of Xe{z) from its standard evolution. 

We thank U. Seljak and M. Zaldariagga for the use of CMBFAST |l|]. R.J.S. was supported in part by the DOE 
(DE-FG02-91ER40690). 
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